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1.0  INTRODUCTION 


1.1  Purpose 

Z-fiber  reinforcement  is  being  examined  for  its  ability  to  increase  the  through-the- 
thickness  strength  of  advanced  composite  laminates,  joints  and  structural  interfaces.  If  this 
technology  proves  successful,  designers  will  be  able  to  use  it  to  help  solve  problems 
associated  with  cocured  composite  structure  including  skin-stiffener  separation  [1  &  2] 
and  hat  stiffener  runout  [3],  Also,  it  would  find  application  in  the  development  of  damage 
tolerant  skins  [4  &  5],  survivable  composite  wing  structure  [6]  and  unitized  aircraft 
structure  [7],  Obtaining  a  fundamental  understanding  of  the  interaction  between  the  z-fiber 
reinforcement  and  the  structure  in  which  it  is  embedded  is  key  to  the  development  of  this 
technology.  To  study  this  interaction,  micromechanical  stress  analyses  were  performed 
based  on  the  finite  element  method  (FEM).  The  results  of  these  analyses  are  reported 
herein. 

1.2  Background 

Laminated  fiber  reinforced  materials  have  high  in-plane  strength  to  weight  ratios. 
Because  of  this,  they  are  the  materials  of  choice  for  the  construction  of  the  wings  and 
fuselage  of  fighter  aircraft.  A  preferred  form  of  laminated  construction  is  integrally 
stiffened  cocured  structure  because  this  method  of  fabrication  reduces  fabrication  and 
assembly  costs.  The  drawback  to  the  use  of  laminates  and  of  cocured  structure  is  that 
laminates  have  poor  through-the-thickness  strength.  Moderate  out-of-plane  loads  can  lead 
to  delaminations,  skin-stiffener  separations  and  failed  joints. 

In  an  attempt  to  improve  the  through-the-thickness  properties  of  laminates, 
researchers  have  examined  the  use  of  stitching,  toughened  resin  systems  and  z-fiber 
reinforcement.  Stitching  involves  the  sewing  of  fibers  through  the  thickness  of  laminates 
and  through  structural  interfaces.  The  stitches  directly  resist  pull-off  loads  and  joint 
reactions.  The  drawbacks  to  stitching  include  a  substantial  capital  investment,  a 
degradation  of  laminate  in-plane  strength  and  stiffness  properties,  and  labor  intensive 
construction  procedures  for  large  structural  parts.  In  toughened  resin  systems,  the 
formulation  of  the  resin  is  changed  in  order  to  improve  the  through-the-thickness 
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structural  properties.  Toughened  resin  systems,  however,  are  costly  and  still  do  not 
provide  adequate  fracture  toughness.  Z-fiber  reinforcement  involves  the  direct  insertion  of 
reinforcing  fibers  in  the  through-the-thickness  direction  of  laminates  and  structural 
interfaces.  (In  analytical  models,  the  through-the-thickness  coordinate  line  is  usually 
designated  the  “z-axis”.  The  term  “z-fibers”  derives  from  the  orientation  of  the  z-fibers 
along  this  axis.)  At  this  point  in  time,  z-fiber  reinforcement  appears  to  offer  many  benefits. 

The  z-fiber  reinforcement  process  was  developed  by  Foster-Miller,  Inc.  of  Waltham, 
Massachusetts.  Their  efforts  were  funded  by  in-house  investment  and  through  various 
Small  Business  Innovative  Research  programs  (SBIR).  Recently,  Foster-Miller  announced 
the  formation  of  Aztex,  Inc.;  a  majority  owned  subsidiary  of  the  parent  company.  Aztex 
will  serve  as  the  supplier  of  z-fiber  products,  equipment  and  engineering  services. 

The  typical  z-fiber  reinforcing  process  [8]  starts  by  placing  a  release  film,  a  z-fiber 
preform  and  a  rigid  tool  onto  a  laid-up  prepreg  laminate.  The  z-fiber  preform  consists  of  a 
structural  foam  that  contains  the  reinforcing  fibers.  The  part  is  then  vacuum  bagged  and 
processed  using  a  standard  autoclave  cure  cycle.  The  heat  softens  the  preform  which 
collapses  under  the  applied  pressure.  The  z-fiber  reinforcement  is  thereby  driven  into  the 
laminate.  (The  foam  is  chosen  for  its  ability  to  provide  lateral  support  to  the  fibers  during 
the  insertion  process.)  Upon  removal  from  the  autoclave,  the  compacted  foam  is  removed 
and  discarded.  The  process  is  completed  with  the  removal  of  any  pin  material  that  projects 
above  the  laminate. 

Z-fiber  preforms  are  made  with  a  specially  designed  insertion  machine  that  embeds  the 
reinforcing  fibers  into  a  structural  foam.  The  insertion  machine  can  be  programmed  to 
insert  the  z-fibers  to  various  lengths,  spacing  and  patterns  [9],  Preforms  are  available  in  a 
variety  of  thicknesses  that  range  from  1  mm  to  19  mm  for  single  fiber  insertions  and  up  to 
51  mm  for  multiple  fiber  insertions.  The  maximum  areal  densities  are  1%  for  standard 
tape,  3%  for  standard  fabric  and  7.5%  for  modified  fabric.  A  number  of  different 
reinforcing  materials  are  in  use  including  pan  and  pitched  based  carbon/epoxy,  carbon 
bismaleimides,  glass/epoxy,  titanium,  refractory  alloys  and  stainless  steel.  The  diameters  of 
the  reinforcing  fibers  range  from  0. 15  mm  to  1  mm. 
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The  types  of  prepreg  systems  to  which  z-fiber  reinforcement  is  applied  include 
AS4/3501-6,  IM7/8551-7,  PEEK/APC-2,  and  S-2(glass)/polyester.  Consolidated 
thermoplastics  have  also  been  reinforced  with  an  ultrasonic  insertion  method. 

Foster-Miller  claims  that  the  z-fibers  benignly  penetrate  laminations  when  proper 
processing  procedures  are  employed.  There  is  minimal  fiber  breakage.  This  claim  is 
supported  by  photomicrographs  that  show  that  the  ply  fibers,  that  were  originally  in  the 
path  of  the  penetrating  z-fibers,  gently  bend  around  the  z-fibers.  This  displacement  of  the 
fibers  results  in  the  formation  of  elliptically  shaped  resin  pools  (see  Figure  1).  On  the 
minor  axis  of  these  ellipses,  the  ply  fibers  are  in  close  proximity  to  the  z-fibers.  Along  the 
major  axis,  the  resin  pool  can  extend  for  several  z-fiber  radii.  The  orientations  of  the 
ellipses  change  from  ply  to  ply  as  the  individual  ply  orientations  change. 

Test  results  [9]  indicate  that  z-fiber  reinforcement  leads  to  an  18  fold  increase  in 
Mode  I  fracture  toughness.  (The  definition  of  this  material  property  does  not  strictly  apply 
to  crack  growth  in  z-fiber  reinforced  laminates.  The  test  results  must  be  interpreted  on  a 
structural  level  as  opposed  to  a  material  level.)  The  use  of  z-fiber  reinforcement  also  leads 
to  a  50%  increase  in  compression  after  impact  strength  over  unreinforced  laminates. 
Reinforced  laminates  do  not  experience  any  decrease  in  compression  strength  and  retain 
from  91  to  98%  of  their  unreinforced  tensile  strength  and  modulus.  A  50%  improvement 
in  ultimate  stiffener  pull-off  loads  has  been  recorded. 

A  substantial  SBIR  investment  has  been  made  in  z-fiber  technology.  Although  some 
of  this  money  was  used  to  support  process  development,  a  large  portion  of  it  was  used  to 
develop  and  test  various  structural  configurations.  Currently,  two  large  scale  technology 
development  programs  are  making  use  of  z-fiber  reinforcement;  Lockheed’s  ccRobust 
Composite  Sandwich  Structures”  program  and  McDonnell-Douglas’s  “Advanced 
Lightweight  Aircraft  Fuselage  Structures”  program  [7],  The  emphasis  of  these  programs  is 
on  the  demonstration  of  new  structural  technologies.  These  two  programs  and  the  SBIR 
efforts  do  not  include  tasks  to  investigate  the  fundamental  mechanical  behavior  of  z-fibers. 
A  fiiture  Great  Lakes  Composite  Consortium  (GLCC)  program,  directed  specifically  at  the 
development  of  z-fiber  reinforcement  technology,  will  include  an  emphasis  on  engineering 
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science.  However,  early  feed  back  from  the  prime  airframers  shows  that  they  will  be  vying 
for  making  and  testing  the  structural  parts  that  give  them  their  greatest  concerns. 

What  is  missing  in  these  research  and  development  activities  is  a  basic  investigation  of 
how  the  z-fibers  interact  with  the  structure  in  which  they  are  embedded.  The  success  of 
this  technology  is  dependent  on  a  clear  understanding  of  the  mechanics  of  this  interaction. 

In  a  typical  aircraft  stress  analysis  that  is  based  on  the  finite  element  method,  the 
resolution  of  the  models  is  rarely  finer  than  the  skin  thicknesses  or  the  individual  ply 
thicknesses.  The  mechanics  of  z-fiber  reinforcement  has  been  examined  to  this  level  by 
Foster-Miller.  Proposals  to  the  GLCC  program  are  expected  to  include  the  same  kind  of 
modeling  with  individual  fibers  represented  by  single  finite  elements.  At  this  point,  the 
suspicion  is  that  these  models  will  not  be  able  to  fully  characterize  the  interaction  between 
the  z-fiber  and  the  laminate.  If  the  stress  transfer  between  the  z-fibers  and  the  surrounding 
structure  follows  a  shear  lag  type  of  distribution,  then  highly  local  stress  concentrations 
will  develop  around  the  z-fibers  at  any  free  boundary.  The  characteristic  dimension  of 
these  stress  concentrations  is  expected  to  be  less  than  a  typical  ply  thickness.  Therefore 
the  proposed  level  of  finite  element  modeling  will  not  be  able  to  resolve  the  magnitude  and 
distribution  of  these  stress  concentrations.  However,  it  is  this  type  of  information  that  will 
be  critical  to  the  development  of  the  optimum  form  and  usage  of  z-fiber  reinforcement. 

To  address  this  issue,  the  Strength  Branch  of  the  Naval  Air  Warfare  Center  -  Aircraft 
Division,  Warminster  (NAWCADWAR),  under  the  sponsorship  of  the  Office  of  Naval 
Research,  undertook  an  effort  to  examine  the  stress  fields  that  are  generated  in  laminates 
by  the  introduction  of  z-fibers.  The  approach  was  based  on  the  development  of  finite 
element  models  that  have  a  fine  enough  resolution  to  examine  stress  concentrations  about 
the  z-fibers.  In  addition,  less  refined  models  were  used  to  examine  the  effect  of  z-fibers  on 
crack  tip  stress  fields. 
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2.0  ANALYSIS 


2.1  Problem  Descriptions 

Two  types  of  structural  problems  were  addressed.  The  first  involved  the  application 
of  thermal  loads  to  z-fiber  reinforced  plates.  The  second  examined  the  effect  of  z-fibers  on 
crack  initiation  and  growth. 

Thermal  stresses  are  generated  when  materials  with  different  coefficients  of  thermal 
expansion  (COTE)  are  constrained  to  respond  in  unison  to  thermal  loads.  In  general,  the 
COTE  of  z-fibers  and  of  the  laminates  in  which  they  are  embedded  are  different.  Thus, 
both  in-plane  and  through-the-thickness  stresses  are  generated  when  z-fiber  reinforced 
laminates  experience  a  change  in  temperature.  The  magnitude  and  distribution  of  these 
stresses  are  of  concern  for  thermal  fatigue  and  for  overloads  resulting  from  thermal  loads 
acting  in  conjunction  with  mechanical  loads. 

Z-fiber  reinforcement  can  also  affect  the  growth  and  development  of  in-plane  cracks. 
For  crack  initiation,  the  z-fibers  can  help  to  shield  the  laminate  interfaces  from  damaging 
stresses.  Also,  when  in-plane  cracks  have  penetrated  throughout  a  region  that  contains  an 
intact  z-fiber,  the  z-fiber  can  transmit  loads  between  the  separated  crack  surfaces.  This  will 
reduce  the  stresses  at  the  tips  of  the  crack  and  will  arrest  the  crack  growth.  To  produce 
these  benefits,  the  z-fibers  receive  load  from  the  laminate  through  interfacial  shear 
stresses.  These  stresses,  which  will  be  concentrated  around  the  z-fibers  in  the  vicinity  of 
the  separated  surfaces,  can  introduce  new  failure  modes  (e.g.,  through-the-thickness 
cracks,  fiber  slippage,  surface  cracks,  etc  ). 

To  examine  the  thermal  load  and  the  crack  stress  issues,  a  z-fiber  reinforced  plate  was 
analyzed.  The  plate  consisted  of  an  18  ply  quasi-isotropic  lamination  of  AS4/3501-6.  The 
z-fibers  were  taken  to  be  T300/5208  circular  pins  with  a  diameter  of  .280  mm.  Table  1 
lists  the  mechanical  properties  of  these  materials  as  well  as  those  of  pure  resin  and 
titanium.  The  total  plate  thickness  was  3.50  mm.  The  areal  density  of  the  z-fibers  was  .5%, 
with  the  z-pins  arranged  in  a  square  array  with  a  spacing  of  3.50  mm. 

The  z-fiber  reinforced  plate  was  analyzed  using  the  axisymmetric  and  the  three- 
dimensional  stress  models  shown  in  Figure  2. 
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2.2  Axisymmetric  Stress  Model 

In  the  axisymmetric  model,  a  single  z-fiber,  the  surrounding  resin  pool  and  the 
laminate  were  modeled  as  three  concentric  cylinders.  The  radius  of  the  model  was  1.75 
mm  which  is  half  the  distance  to  the  adjacent  z-fibers.  The  resin  pool  was  assumed  to  have 
a  thickness  of  .030  mm. 

Two  assumptions  are  implicit  in  this  model.  The  first  is  that  the  individual  ply 
orientations  are  not  critical  to  the  stresses  of  concern  (the  through-the-thickness  normal 
and  shear  stresses).  Based  on  this,  the  laminate  was  treated  as  a  monolithic  quasi-isotropic 
material.  The  second  assumption  is  that  the  interfacial  stresses  between  the  z-fiber  and  the 
surrounding  resin  pool  are  relatively  insensitive  to  the  pool  thickness.  Therefore  a  circular 
pool,  uniform  throughout  the  thickness,  was  used  to  represent  an  elliptically  shaped  pool 
whose  major  axis  changes  orientation  with  each  change  in  ply  direction.  (Subsequent  work 
will  be  needed  to  validate  this  last  assumption.  The  work  reported  here  is  a  reasonable 
attempt  to  determine  the  mechanical  response.) 

The  axisymmetric  model  was  used  to  examine  the  thermal  stresses  that  are  induced  by 
a  uniform  temperature  change  in  the  plate.  Also,  it  was  used  to  determine  the 
axisymmetric  crack  tip  stresses  that  form  in  a  region  about  a  z-fiber  under  a  uniform 
through-the-thickness  crack  opening  load  (see  Figure  3).  Because  the  model  and  the 
applied  loads  were  axisymmetric,  the  model  was  analyzed  with  a  two-dimensional 
axisymmetric  finite  element  mesh  (see  Figure  4).  The  mid-plane  symmetry  of  the  plate  was 
used  to  reduce  the  size  of  this  mesh.  (The  mesh  extends  from  the  mid-depth  to  the  surface 
of  the  plate.)  Figure  4  shows  the  X-coordinate  line  extending  from  the  centerline  of  the  z- 
fiber  into  the  radial  direction.  The  Y-coordinate  extends  from  the  mid-depth  of  the  plate 
into  the  axial  direction  of  the  z-fiber.  The  notation  on  this  figure  indicates  the  location  of 
the  z-fiber  (MP=1),  the  resin  (MP=2)  and  laminate  (MP=3).  The  mesh  density  was  sized 
to  analyze  the  stress  concentration  that  would  occur  in  the  interfacial  area  between  the  z- 
fiber  and  the  resin  when  thermal  loads  are  applied.  An  enlarged  view  of  the  interfacial  area 
at  the  surface  of  the  model  is  shown  in  Figure  5.  The  aspect  ratios  of  the  elements  have  an 
upper  limit  of  approximately  3  to  1 . 
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2.2.1  Thermal  Stress  Analysis 

Figures  6  and  7  show  the  displacements  that  occur  in  the  axisymmetric  finite  element 
mesh  for  a  10  °C  rise  in  temperature.  Numerical  values  for  the  radial  (Disp  X)  and  axial 
(Disp  Y)  displacements  are  shown  in  Figures  8  and  9.  The  displacements  in  the  axial 
direction  on  the  outer  edge  of  the  model  match  the  displacements  that  would  be  found  for 
the  general  laminate.  This  establishes  the  validity  of  using  the  model  to  analyze  the  thermal 
behavior  of  the  plate.  The  displacements  of  the  z-fiber  at  the  surface  of  the  model  are 
approximately  the  value  predicted  by  a  Rule-of-Mixtures  model.  The  radial  displacements 
in  the  z-fiber  are  seen  to  be  negative.  This  is  caused  by  the  particular  combination  of 
materials  used  in  this  model.  The  laminate,  which  has  a  high  in-plane  modulus  and  a  low 
COTE,  acts  as  a  barrier  to  the  thermal  radial  growth  of  the  z-fiber  and  the  resin.  The 
combination  of  the  moduli  and  the  COTE  of  the  resin  dominates  the  corresponding  in¬ 
plane  properties  of  the  z-fiber.  The  result  is  that  the  resin  pushes  inwards  from  the 
laminate  into  the  z-fiber.  The  displacement  plots  also  show  the  intense  shearing  that 
occurs  near  the  model’s  upper  surface  in  the  elements  that  are  used  to  model  the  resin 
pool. 

The  through-the-thickness  normal  stresses  (Sigma  Y)  are  shown  in  Figures  10  and 
1 1 .  The  stresses  in  the  z-fiber  are  tensile  since  the  laminate’s  greater  thermal  growth  pulls 
on  the  z-fiber.  On  the  other  hand,  the  stresses  in  the  laminate  are  compressive  since  the  z- 
fiber  retards  the  laminate’s  thermal  movements.  The  laminate’s  compressive  stresses  tend 
to  zero  along  the  radial  coordinate. 

The  radial  stresses  (Sigma  X)  are  shown  in  Figures  12  and  13.  As  the  radial 
displacement  field  indicates,  the  radial  stresses  in  the  z-fiber,  the  resin  pool  and  the 
adjacent  laminate  are  compressive.  The  radial  stresses  in  the  laminate  tend  to  zero  in  the 
radial  direction.  The  hoop  stresses  (Sigma  Z)  are  shown  in  Figures  14  and  15.  These 
stresses  are  compressive  in  the  z-fiber  and  in  the  resin  pool  but  become  tensile  in  the 
adjacent  laminate.  The  hoop  stress  also  tends  to  zero  in  the  radial  direction.  Figures  13 
and  1 5  show  how  the  stresses  deviate  from  these  patterns  because  of  the  free  surface. 
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The  through-the-thickness  shear  stresses  are  shown  in  Figures  16  and  17.  The  highly 
stressed  area  is  seen  to  occur  in  the  interface  between  the  z-fiber  and  the  resin  pool  just 
below  the  free  boundary. 

2.2.2  Axisymmetric  Crack  Tip  Stress  Analysis 

A  crack  was  introduced  at  mid-depth  in  the  axisymmetric  model.  The  crack 
penetrated  from  the  outer  rim  of  the  model  inwards  along  the  radial  direction  (see  Figure 
3).  The  crack  was  introduced  into  the  finite  element  mesh  by  releasing  the  restraints 
against  the  Y-direction  displacements  along  the  base  of  the  mesh.  Figure  18  shows  the 
mesh  and  the  boundary  conditions  for  a  crack  that  has  penetrated  to  a  radial  location  of 
.700  mm,  that  is,  to  within  2  z-fiber  diameters  of  the  face  of  the  z-fiber. 

Figure  19  shows  the  displacements  that  occur  in  the  mesh  for  crack  opening  loads 
applied  around  the  rim  of  the  model  (18  nodal  forces  of  .00  IN  acting  in  the  vertical 
direction  (see  Figure  18)).  Figures  20  through  23  show  the  associated  stress  fields. 

As  a  basis  for  comparison,  the  axisymmetric  model  was  used  to  analyze  a  section  of 
laminate  that  contains  no  z-fiber  reinforcement.  For  this  case,  the  parts  of  the  mesh,  that 
were  previously  used  to  represent  the  z-fiber  and  the  resin  pool,  were  reassigned  the 
material  properties  of  the  laminate.  The  crack  opening  loads  and  the  displacement 
boundary  conditions  remained  the  same.  The  resulting  displacement  field  is  shown  in 
Figure  24.  The  through-the-thickness  normal  stresses  are  shown  in  Figure  25.  The  stress 
components  in  the  element  at  the  tip  of  the  crack  for  the  reinforced  and  unreinforced  cases 
are  listed  and  compared  in  Table  2.  A  comparison  between  the  through-the-thickness 
stresses  that  occur  along  the  mid-plane  of  the  reinforced  (Series  1)  and  the  unreinforced 
(Series  2)  cases  is  shown  in  Figure  26. 

The  final  analysis  performed  on  the  axisymmetric  model  considered  a  crack  that  had 
penetrated  up  to  the  face  of  the  z-fiber.  The  crack  opening  loads  remained  the  same.  The 
resulting  displacements  are  shown  in  Figure  27.  The  stress  fields  are  shown  in  Figures  28 
and  29. 

2.3  Three  Dimensional  Stress  Model 

Mode  I  crack  opening  loads  were  applied  along  the  entire  edge  of  the  plate  from 
which  the  three-dimensional  stress  model  was  taken.  A  mid-plane  crack  on  this  edge  was 
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assumed  to  extend  inwards  to  a  point  midway  between  the  first  and  second  z-fiber  (see 
Figure  30).  Given  the  scale  of  the  z-fiber  spacing,  the  mechanical  response  to  the 
described  load  in  each  row  of  z-fibers  will  be  approximately  the  same.  (The  only  exception 
would  be  for  rows  near  the  unloaded  lateral  edges,  the  analysis  of  which  is  not  addressed 
here). 

To  analyze  the  stresses  in  the  crack  tip,  a  three-dimensional  stress  model  was 
formulated  which  consisted  of  a  rectangular  block  of  laminate  containing  two  z-fibers  (see 
Figures  2  and  30).  The  block  was  taken  from  the  loaded  end  of  the  plate  with  the  long 
edge  of  the  block  directed  towards  the  plate’s  interior.  The  width  of  the  block  was  3.50 
mm  and  its’  length  was  7.00  mm.  The  z-fibers  were  located  along  the  widthwise  centerline 
of  the  block  with  the  z-fibers  centered  1.75  mm  from  the  front  and  back  faces  of  the  block. 
As  defined,  the  model  has  two  planes  of  structural  symmetry.  Since  the  crack  tip  stresses 
were  of  interest,  the  pools  of  resin  about  the  z-fibers  were  not  included  in  the  model. 

To  achieve  the  quasi-plane  strain  conditions  that  exist  between  the  rows  of  z-fibers, 
the  out-of-plane  displacements  on  the  side  faces  of  the  model  were  held  fixed.  At  the  rear 
of  the  model,  the  displacements  were  assumed  to  be  free  of  any  restraint.  Implicit  in  this 
assumption  is  that  the  stiffness  of  the  plate  past  the  second  z-fiber  has  a  negligible  affect 
on  the  crack  tip  stresses. 

The  three-dimensional  model  was  analyzed  with  a  mesh  of  8  noded  isoparametric 
solid  elements.  Because  the  structure  and  the  loading  were  symmetric  about  the  mid-plane 
of  the  plate  and  about  a  vertical  plane  through  the  centers  of  the  two  z-fibers,  only  the 
upper  right  hand  quarter  of  the  stress  model  was  analyzed.  Figure  3 1  shows  the  mesh  with 
the  origin  of  the  reference  coordinate  system  located  at  the  point  where  the  two  planes  of 
structural  symmetry  intersect  the  plane  perpendicular  to  the  tip  of  the  crack.  The  X- 
coordinate  line  of  this  system  extends  towards  the  free  edge  of  the  plate,  the  Z-coordinate 
line  extends  upwards  through  the  thickness  of  the  plate,  and  the  Y  coordinate  line 
completes  a  right  handed  coordinate  system.  The  notation  on  Figure  31  indicates  the 
location  of  the  z-fiber  (MP=1)  and  the  laminate  (MP=2). 

The  intention  of  the  analysis  was  to  examine  the  effect  that  z-fibers  have  on  crack 
propagation.  This  was  accomplish  by  comparing  a  case  in  which  the  z-fiber  remains  intact 
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in  the  cracked  region  against  a  case  were  the  z-fiber  is  severed  in  the  plane  of  the  crack. 
Mechanically,  the  last  case  is  equivalent  to  no  reinforcement  being  present  at  all.  The 
boundary  conditions  used  in  the  first  case  are  shown  in  Figures  32  and  33.  Note  that  the 
intact  z-fiber  was  modeled  by  fixing  the  displacements  in  the  plane  of  the  crack.  The  right 
hand  side  of  Figure  32  shows  the  applied  pressure  (1  MPa).  The  boundary  conditions  in 
the  second  case  were  identical  to  the  first  with  the  exception  that  the  z-fiber  displacements 
were  released  at  the  crack  (see  Figure  34).  For  making  this  comparison,  the  exact  level  of 
stress  concentration  about  the  tip  of  the  crack  was  not  needed  so  that  the  three 
dimensional  mesh  used  here  is  not  as  refined  as  the  axisymmetric  mesh. 

2.3.1  Crack  Tip  Stress  Analysis 

Figure  35  shows  the  displacements  that  occur  when  the  crack  opening  loads  are 
applied  and  the  z-fiber  remains  intact.  Figures  36  and  37  show  different  views  of  this 
displacement  field.  The  associated  bending  stresses  (Sigma  X),  the  through-the  thickness 
normal  stresses  (Sigma  Y),  the  in-plane  lateral  stresses  (Sigma  Z)  and  the  shear  stresses 
(Tau_XZ)  are  shown  in  Figures  38  through  41.  The  tractions  at  the  rear  of  the  model  are 
negligible,  which  validates  the  assumption  to  neglect  the  stiffness  of  the  adjacent  laminate. 
The  displacements  and  stresses  for  when  the  z-fiber  was  severed  are  shown  in  Figures  42 
through  47.  A  comparison  between  the  two  cases  for  the  peak  stresses  at  the  tip  of  the 
crack  is  made  in  Table  3. 
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3.0  CLOSURE 


3.1  Discussion  and  Conclusions 

In  the  thermal  stress  analysis  it  was  shown  that  thermal  loads  lead  to  stress 
concentrations  in  an  AS4/3501-6  laminate  reinforced  with  T-300  z-fibers.  The  maximum 
shear  stress  in  the  interfacial  resin  was  shown  to  be  3.52  MPa  for  a  temperature  rise  of  10 
°C  (see  Figure  16).  Since  an  aircraft  component  can  routinely  experience  temperature 
changes  of  50  to  60  degrees  and  since  a  typical  resin  has  a  shear  strength  of  60  MPa,  the 
predicted  thermal  stress  concentrations  are  not  negligible. 

Matching  the  COTE  of  the  z-fibers  to  that  of  the  laminate  will  help  to  control  the 
thermally  induced  shear  stresses  For  instance,  titanium  would  be  a  better  candidate  for 
reinforcing  an  AS4/3501-6  laminate  since  its  COTE  is  higher  than  that  of  T-300.  Note 
however,  that  this  higher  COTE  and  titanium’s  large  elastic  modulus  may  lead  to  the 
build-up  of  in-plane  stresses.  Therefore,  unless  the  materials  of  construction  are  perfectly 
matched,  thermal  stresses,  thermal  cycling  and  thermal  fatigue  are  of  concern.  These 
issues  should  be  addressed  in  the  environmental  test  plans  of  all  of  those  who  are 
evaluating  the  use  of  this  technology. 

The  stress  analyses  showed  that  the  z-fibers  resist  through-the-thickness  loads  only 
when  those  loads  are  applied  in  the  immediate  vicinity  of  a  z-fiber  (see  Figure  26  and 
Table  2).  This  means  that,  for  typical  z-fiber  areal  densities,  most  of  a  laminate  will  remain 
unprotected  from  through-the-thickness  loads.  Therefore  z-fibers  can  do  very  little  to 
prevent  crack  initiation. 

Once  an  in-plane  crack  has  propagated  throughout  a  region  that  contains  an  intact  z- 
fiber,  the  z-fiber  will  pick-up  the  through-the-thickness  loads.  This  will  greatly  reduce  the 
stresses  at  the  tip  of  the  crack  (see  Table  3  and  compare  Figures  39  and  45).  Therefore  z- 
fibers  have  a  significant  affect  on  crack  propagation. 

When  z-fibers  resist  through-the-thickness  stresses,  there  is  a  shear  lag  type  of  stress 
transfer  between  the  z-fiber  and  the  surrounding  material  (see  Figure  29).  The  magnitude 
of  this  stress  concentration  needs  to  be  determined  in  order  to  avoid  unexpected  failure 
modes  and  to  achieve  an  optimized  design 
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The  characteristic  dimension  of  the  stress  concentrations  examined  in  this  work  are 
less  than  the  individual  composite  ply  thicknesses.  Aircraft  analysts  rarely  develop  stress 
models  as  fine  as  ply  thicknesses.  Therefore  current  modeling  methods  will  not  be  able  to 
predict  the  local  effects  of  z-fiber  usage  Also,  in  order  to  avoid  computational  expense, 
there  is  a  need  to  develop  a  micromechanical  model  that  can  bridge  the  gap  between  the 
typical  aircraft  stress  analysis  and  the  refined  FEM  meshes  presented  in  this  work. 
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T300/5208  Z-Fiber  (Transversely  Isotropic) 


Axial  Modulus 
Transverse  Modulus 
Axial  Poisson’s  Ratio 
Axial  Shear  Modulus 
Transverse  Shear  Modulus 
Axial  COTE 
Transverse  COTE 


152.  GPa 
10.3  GPa 
0.31 

5.70  GPa 
3.60  GPa 
0.54  ps/°C 
19.8  pe/°C 


Resin  (Isotropic) 

Extensional  Modulus 
Poisson’s  Ratio 
Shear  Modulus 
COTE 

Tensile  Strength 
Compressive  Strength 
Shear  Strength 


3.52  GPa 
0.38 

1.28  GPa 
67.9  p£/°C 
40  MPa 
69  MPa 
60  MPa 


AS4/3501-6  Laminate  (Quasi-Isotropic) 


In-Plane  Modulus 
Through-the-Thickness  Modulus 
Through-the-Thickness  Poisson’s  Ratio 
In-Plane  Shear  Modulus 
Through-the-Thickness  Shear  Modulus 
In-Plane  COTE 
Through-the  Thickness  COTE 


54.3  GPa 

10.6  GPa 
.0508 
3.45  GPa 

20.7  GPa 
1.93pe/°C 
32.0  ps/°C 


Titanium  (Isotropic) 

Extensional  Modulus 
Poisson’s  Ratio 
Shear  Modulus 
COTE 


110.  GPa 
0.30 

4.34  GPa 
8.82  (J.e/°C 


Table  #1  Material  Properties 
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Stress  Component  Severed  Z-Fiber  Intact  Z-Fiber  Reduction  of  Stress 
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Table  #3  Reduction  in  Nodal  Stress  due  to  Intact  Z-Fiber 


Photomicrographs  of  Z-Fiber  Reinforced  Laminate  (160x) 
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Quasi-Isotropic  Laminate 
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Figure  6  — Thermally  Induced  Displacements 
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Figure  7  —  Enlarged  View  of  Thermally  Induced  Displacements 
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Figure  9  —  Axial  Displacements  (M) 
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Figure  10  —  Through-the-Thickness  Normal  Stress  Field  (Pa),  Thermal  Loads 
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Figure  11  —  Enlarged  View  of  the  Through-the-Thickness  Normal  Stress  Field  (Pa) 
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Figure  12  —  Radial  Stress  Field  (Pa),  Thermal  Loads 
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Figure  13  —  Enlarged  View  of  the  Radial  Stress  Field  (Pa)?  Thermal  Loads 


Figure  14  —  Hoop  Stress  Field,  (Pa),  Thermal  Loads 
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Figure  15  —  Enlarged  View  of  the  Hoop  Stress  Field  (Pa),  Thermal  Loads 
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Figure  17  —  Enlarged  View  of  the  Shear  Stress  Field  (Pa),  Thermal  Loads 
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Figure  19  —  Displacements  due  to  Crack  Opening  Loads 
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Figure  20  —  Through-the-Thickness  Normal  Stress  Field  (Pa),  Crack  Opening  Loads 


-49QE+Q05 


NAWCADWAR-95-1 2-TR 


Lnininmin^^r^^^ 
©00B000O00 
e>©©©©©©©©© 
+  +  +  +  -!-  +  +  +•+  + 

LiJ  LiJ  LlJ  LU  LjJ  LU  LlJ  LlJ  LU  LlI 

©©e>©0©©©^cn 

wMro^tin^-cor'BOo 

(\|  (J)  (D  M  0  U)  h  00  •  ■ 

. -  -  •  *  *  n 


<o 

co 

LU 

C£ 


CO 


37 


Figure  21  —  Radial  Stress  Field  (Pa),  Crack  Opening  Loads 


Figure  22  —  Hoop  Stress  Field  (Pa),  Crack  Opening  Loads 
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Figure  23  —  Shear  Stress  Field  (Pa),  Crack  Opening  Loads 
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Figure  25  — Through-the-Thickness  Normal  Stress  Field  (Pa),  Unreinforced  Case 
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Figure  28  —  Through-the-Thickness  Normal  Stress  Field  (Pa),  Crack  Opening  Loads 
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Figure  29  —  Shear  Stress  Field  (Pa),  Crack  Opening  Loads 


Figure  30  —  Mode  I  Loads  and  Edge  Crack  (Not  to  Scale) 
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Figure  31  —  Three-Dimensional  Finite  Element  Mesh 
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Figure  32  —  Boundary  Conditions  as  Viewed  from  the  Negative  Y-Axis,  Pin  Intact 
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Figure  33  —  Boundary  Conditions  as  Viewed  from  the  Positive  Z-Axis,  Pin  Intact 
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Figure  34  —  Boundary  Conditions  as  Viewed  from  the  Negative  Y-Axis,  Pin  Severed 
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Figure  35  —  Displacements  due  to  Crack  Opening  Loads,  Pin  Intact 
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Figure  37  —  Displacements  due  to  Crack  Opening  Loads,  Side  View 
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Figure  38  —  Bending  Stresses  (Pa),  Pin  Intact 
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Figure  39  —  Through-the-Thickness  Stresses  (Pa),  Pin  Intact 
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Figure  40  —  In-Plane  Lateral  Stresses  (Pa),  Pin  Intact 
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Figure  41  —  Shear  Stresses  (Pa),  Pin  Intact 
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Figure  42  —  Displacements  due  to  Crack  Opening  Loads,  Pin  Severed 
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Figure  43  —  Displacements  due  to  Crack  Opening  Loads,  Side  View 


Figure  44  —  Bending  Stresses  (Pa),  Pin  Severed 
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